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ON NON-PRIMORDIAL DEUTERIUM PRODUCTION BY ACCELERATED PARTICLES 
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ABSTRACT 

Deuterium plays a crucial role in cosmology because the primordial D/H abundance, in the context 
of big bang nucleosynthesis (BBN) theory, yields a precise measure of the cosmic baryon content. 
Observations of D/H can limit or measure the true primordial abundance because D is thought to 
be destroyed by stars and thus D/H monotonically decreases after BBN. Recently, however, Mullan 
& Linsky have pointed out that D arises as a secondary product of neutrons in stellar flares which 
then capture on protons via n + p — > d + 7, and that this could dominate over direct D production 
in flares. Mullan & Linsky note that if this process is sufficiently vigorous in flaring dwarf stars, it 
could lead to significant non-BBN D production. We have considered the production of D in stellar 
flares, both directly and by n capture. We find that for reasonable flare spectra, n/d < 10 and 
(n + d) / 6 Li < 400, both of which indicate that the n-capture channel does not allow for Galactic D 
production at a level which will reverse the monotonic decline of D. We also calculate the 2.22 MeV 
7-ray line production associated with n capture, and find that existing COMPTEL limits also rule 
out significant D production in the Galaxy today. Thus, we find flares in particular, and neutron 
captures in general, are not an important Galactic source of D. On the other hand, we cannot exclude 
that flare production might contribute to recent FUSE observations of large variations in the local 
interstellar D/H abundance; we do, however, offer important constraints on this possibility. Finally, 
since flare stars should inevitably produce some n-capture events, a search for diffuse 2.22 MeV 7-rays 
by INTEGRAL can further constrain (or measure!) Galactic deuterium production via n-capture. 

Subject headings: nuclear reactions, nucleosynthesis, abundances — cosmology: theory — stars:Hare 
— gamma rays:theory 



1. INTRODUCTION 

In the past two decades big-bang nucleosynthesis 
(BBN) has beco me one of the most important cos- 
mological probes ijOlive. Steigman. fc Walkerll200fi and 
refs. therein). Its success lies in the good agreement 
of predictions for the abundances of four light ele- 
ments D, 3 He, 4 He and 7 Li with their observations. A 
key element is deuterium, which is considered to be 
the best cosmic "baryome ter" among the light elements 
ijSchramm fc Turner! 1 1998j) because of its strong depen- 
dence on the baryon density, or equivalently the baryon- 
to-photon ratio n. Furthermore, now that the cos- 
mic microwave background anisotropy measurements of 
WMAP have independent ly measured n to high pre- 
cision l)Snergel et al.l |2Q03) , D takes a new role as a 
probe of both early universe physics and of astrophysics 
llCvburt. Fields, fc Olivdl2002l l200l . Therefore, it is 
crucial that we fully understand the evolution of deu- 
terium after BBN, in order to correctly infer the pri- 
mordial abundance from observations in the z <C I0 10 
universe. 

Key to D is that there is no significant astrophysical 
production site except for the big bang, and that stars 
destro y D in their fully co nvective, pre-main sequence 
phase (|Bodenheimerl fl9 66) . Together, these guarantee 
that any measurement of D is a solid lower bound on 
the primordial abundance, and that in sufficiently prim- 
itive environments D should be essentially primordial. 
The lack of astrophysic al D production was established 
in the classic paper by lEpstein. Lattimer. fc Schramm! 
l)1976t hereafter ELS), who considered all known sites 



in which any nucleosynthesis occurs, and demonstrated 
that none of these produce D in significant quantities. 

Deuterium is observed in diverse astrophysical settings. 
In high-redshift (z ~ 3) QSO absorptio n systems, neu- 
tral D is observed. The 5 best sy s tems l|Bu rles fc Ty^lerl 
1998a b; Kir kman et all l2003t lO'Meara et alJ I200H 
iPettini fc Bowenl I200H) give (D/H) QSO als = (2.78 ± 
0.29) x 10~ 5 . The D abundance in the solar nebula 
is inferred from solar wind observations of 3 He (which 
measure pre-solar D+ 3 He) min us meteoritic determina- 
tions of 3 He alone. These give l|Geiss & Gloeckler 1998) 
a value (D/H) pre _o = (2.1 ± 0.5) x 10~ 5 , which probes 
proto-solar mater ial 4.6 Gyr ago, or at z ~ 0.4. D/H 
is als o observed ijLinskv et al.l 119931 iRoeerson fc York! 
119731) in the local interstellar medi um (ISM); recent 
FUSE observations l|Moos et al.ll200l give (D/H) IS m = 
(1.52 ± 0.08) x 10" 5 in the Local Bubble today, at z = Q. 

The central ELS argument-that the big bang is the 
sole source of D-has stood the test of time, borne out 
by the drop in D/H from its high-redshift to pre-solar 
to local ISM abundance. But given the crucial role of D 
in cosmology, it is important to carefully examine the 
assumptions made, and to identify any possible loop- 
holes. Moreover, such an effort becomes crucial in light 
of growing evidence that D/H has large variations over 
short distances in the local ISM. Recent FUSE observa- 



tions i Hoooes et al.M2003l) add weight to earlier sugges- 
tions ijVidal-Madiarl 120011 and refs. therein) that D/H 



can vary by as much as a factor of 2 over different lines 
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of sight, pointing to inhomogeneities on scales > 100 pc. 1 
Such a re-analy sis of ELS was recently carried out by 
iMullan fc Linskvl |l999. hereafter ML). Specifically, they 
discussed D production by suprathermal energetic par- 
ticles in flare stars. Although ELS had considered flare 
production, ML note that ELS had neglected flare pro- 
duction of neutrons, which can then go on to produce D 
via radiative capture onto a proton 

n + p — > d + 7 (1) 
with the emission of a photon with energy E~ = 2.223 
MeV. ML do not make a detailed calculation of n yields 
from spallation reactions, but note that n production can 
be large compared to direct d production, depending on 
the energetic particle spectra. ML then suggest that D 
created in such a way can escape from the flare site into 
the stellar wind. Also, because of its low mass D might 
be preferentially ejected relative to the heavier Li, Be, 
an B which would otherwise contaminate the ISM. ML 
thus raise the possibility that flares could be a significant 
source of D in the Galaxy. From the point of view of 
Galactic chemical evolution, this mechanism might then 
explain some or all of the ISM variation in D. Moreover, 
from the point of view of cosmology, the chance of sig- 
nificant D production would call into question the ELS 
argument which underlies the recent spectacular agree- 
ment between high-redshift D/H measurements and the 
predictions of BBN with the WMAP baryon density. 

Thus, the ML scenario has important implications and 
deserves careful further examination. In this paper, we 
expand the analysis of ELS to include the n channel for 
D production suggested by ML. Our approach is to iden- 
tify general nuclear physics constraints that follow from 
detailed calculations of n, d, and 7-ray production by en- 
ergetic particles. The result is to close this loophole-i.e., 
we show that flare D production is not sufficient to offset 
the monotonic decline of D with time. However, flare D 
production should occur at some level, and we find that 
7-ray observations by the recently launched INTEGRAL 
mission may be able to probe this process. 

2. LIGHT ELEMENT PRODUCTION BY FLARES 

Flares are violent releases of magnetic energy during 
which matter reaches temperatures of tens or hundreds 
of million Kelvin, and particles get accelerated to very 
high en ergies. They are associated wit h active stellar 
regions ijHaisch. Strong. fc Rodonol lT991). Nuclear reac- 
tions then occur between flare-accelerated particles and 
the ambient stellar atmosphere. In this section we will 
discuss spallation reactions that result in production of 
neutrons, deuterium and lithium. We will show that, as 
a result of spallation, 1 < n/d < 10 for most spectral 
indices, which means that ML channel is indeed more 
important than channels considered by ELS, confirming 
and quantifying the suggestion of ML. However, we will 
see that even including this channel for d production, 
the concomitant Li production remains large and thus 
severely constrains the possibility that flares are an im- 
portant source of D. 

1 More observations are needed to firmly establish the nature and 
deg ree of D variatio ns, and their (anti-)correlations with metallic- 
ity ISte igman 200ij). Nevertheless, it is already clear that FUSE is 
opening new windows on the Galactic and local evolution of deu- 
terium; this in turn motivates a careful re-examination of the basic 
ELS assumptions which have guided our thinking on these issues. 



2.1. Formalism 

Nuclide production in flares depends on the initial 
spectrum of the flare, its modification as it interacts 
with the stellar atmosphere, and the compositions of 
the flare and the atmosphere. In general, the nucle- 
osynthesis yields in flares have a complicated time de- 
pendence. Fortunately, there are two limiting approxi- 
mations for energetic par ticles in flares: the thin-targe t 
and thick-target models l)Ramatv fe Lingcnfe lteriri975j) . 
In the thin-target model, particle production is assumed 
to occur before the spectrum of accelerated particles is 
modified by ionization energy losses. On the other hand, 
in the thick-target model, energy losses are taken to be 
large, and particles slow due to ionization losses as they 
move downward form the flare region, prior to nuclear 
interactions. 

But are both of these scenarios equally important for 
flare processes? A simple analysis of the mean free path 
of projectile particles sheds light on this question. The 
mean free path of projectiles against nuclear interactions 
is 
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where n is the target number density and a is the reaction 
cross section, while mean free path of a flare particle 
(charge Z, mass number A) that is losing energy due to 
its interaction with the surrounding medium is given by 



A R p {e) 
Z 2 nrrir, 



(3) 



Here, e is the projectile energy per nucleon, n stands for 
the number density of the medium, a is the cross-section 
for a reaction of interest, m p is the proton mass, and R p 
is the ionization energy loss range of protons in units of 
g/cm 2 . In the thin-target limit A <C A c which puts a 
lower limit on the projectile range 
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For low energies, R P (e) ~ 4 x 10~ 4 g/cm 2 (e/MeV) 2 , so 
that eq. J3J implies that the thin-target approximation 
holds for energies 



e > 600 \JZ 2 jA MeV/nucleon 



(5) 



This energy is much larger than those of typical flare 
particles. Thus, we conclude that the thick-target ap- 
proximation is well-satisfied and is the appropriate one 
for flare processes. This result, based on the physics of 
particle propagation, is in agreement with the emp irical 
result of lRamatv. L ingcnfcl ter. fc Kozlovskvl l)1982l) who 
found that flare data are best described by thick-target 
rather than thin-target models. 

In the thick target model, the production of secondary 
par ticles of species I from reaction i + j ' ^ I + ■ ■ ■ is given 
by ijRamatv fc Lingenfelterlll975j) : 



rOO 

Qi = VtVj j 

J 6th 



< dRj 
(m) de 



dc 



N p (e')de' (6) 



where Qi is the total number of secondary particles pro- 
duced per incident flare particle, and eth is the threshold 
energy for a particular reaction. The ambient number 
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density of target-species j is rij and a\j is the cross sec- 
tion for the reaction of interest. Projectile and target 
abundances are given as numbers relative to hydrogen: 
Hi = ni/n p ; these are used to compute the mean target 
mass (m) ~ m p J2i MVil J2j Vj = l-4m p . 

The projectile and target compositions are poorly 
known for flare stars, but one can imagine sev- 
eral possibilities. Since solar flares show en- 
hanced metals and helium relative to hydrogen 
l)Mandzhavidze. Ramatv. fc Kozlovskvlll999|) . one might 
expect similar enrichment for Hare stars. On 
the other hand, one expects M dwarfs in general 

to hav e metallicities similar to those of G or K 

dwarfs dKotone va. Flvnn. Chiappini. fc Matteucclll2002t 
iWvse fc Girmordll995|) . and thus on average to be sub- 
solar in metallicity. Thus we will examine several pos- 
sible compositions and find the impact of composition 
on spallation production. For a fiducial case, we will 
assume solar abundances l) Anders fc Gre vessel 1 1989|) for 
both projectiles an targets. We then will examine the 
effect of more extreme variations: a five-fold increase in 
metals and helium as in solar flares, and a metal-free 
primordial case. 

The spectrum Ni(e) de measures the total number of 
projectile particles with energ y in (e, e + de). Following 
IRamatv fc Lingenfelterl l|1975l 1. we adopt the power-law 
form 

N p {e) = k p e- S (7) 

N i (e)=y i N p (e) (8) 

where k p is the constant determined by normalizing N p 
to 1 proton of energy greater than eo = 30 MeV and s 
is the spectral index which is assumed to be the same 
for all accelerated particles. The choice of this particular 
spectrum normalization is pur ely conventional (see e.g. 
lRamatv"fe Lingenfelter l)1975() h but in any case it will 
not matter in our analysis since we are only interested 
in ratios of productions of different particle species, and 
the normalization will drop out. 

In addition to the particle compositions, the spectral 
index s is the other key parameter for the results we 
present. We will present results for a range of s, but 
some obser vations exist which constrain its value. For 
solar flares. IRamatv. Mandzhavidze. fc Kozlovskvl l)1996f) 
find that 7-ray line ratios in an ensemble of flares point 
to a range of spectral indices of about s = 4 ± 1. 
For flare stars themselves, information on the spectra 
of the ions is not available, but radio data from elec- 
tron synchrotron emission points to electron spectra with 
s e = 2 to 3. Whether the ions should have the same 
indices is, however, unclear. I n the case of the Sun, 
iLin. Mewaldt. fc van Hollebekel l)1982l) find that electron 
spectral indices (for E e < 100 keV) are correlated with 
proton indices, but with s e = 2 corresponding to s ~ 4; 
they find no significant correlation for E e > 200 keV. 
These results hint that the ion spectral indices in flare 
stars may also prefer s ~ 4, but this parameter remains 
uncertain. 

The particle ionization energy loss is measured by 
de_ _ \de/dt\ 

dR pv { ' 

with units of [MeV/(g cm -2 )] where v is the velocity 
of incident particle. We adopt the usual Bethe-Bloch 



Table 1 

Reactions included in our analysis. 



Reaction 






[MeV / nucleon] 


p+p — » n+x 


410.20 


p+ a — > n+x 


22.54 


p+CNO — > n+x 


2.43 


a + a — ► n+x 


9.44 


a +CNO — ► n+x 


0.49 


p+ a — > d+x 


23.07 


p+CNO — > d+x 


17.50 


a + a -» d+ 6 Li 


15.38 


a + C -> d+x 


4.53 


a + O — ► d+x 


5.10 


p+C ^ 6 Li+x 


24.45 


p+N -^ 6 Li+x 


17.52 


p+O ^ 6 Li+x 


23.57 


a + a — > 6 Li +x 


11.19 


o+C -^ 6 Li +x 


7.91 


Q+N -^ 6 Li +x 


2.83 


a+O -^ 6 Li +x 


6.02 



Reference 



Ramatv & Lingenfelter (1975) 
^gmaty^^^jingcnfclter ^^75) 
Ramatv & Lingenfelter (1975) 
Ramatv & Lingenfelter (1975) 
Ramatv & Lingenfelter (1975) 



Meyer (19721 

Meyer (19721 
Mercer et al. (20011 
Bertrand. Peelle. fc Kalbach-Cline (19741 
Bertrand. Peelle. fc Kalbach-Cline (19741 

Read fc Viola (19841 
Read fc Viola (19841 
Read fc Viola (19841 
Mercer et al. (20011 
Read fc Viola (19841 
Read fc Viola (19841 
Read fc Viola (19841 



energy formula (|BethellI930|) 

de_ 47rzZ 2 e 4 
dR A(m)m e v 



In 



2'y 2 m e v z 
I 



(10) 



In this equation Z is the charge of projectile, z is number 
of electrons per atom which is approximately equal to 1 
in the solar atmosphere, A is the number of nucleons of 
projectile, m e is electron mass, and the mean excitation 
potential is the value for hydrogen, / = 13.6 eV. 

Table 1 lists the reactions we have included, their lab- 
frame threshold energies, and references for the cross- 
sections. We note here that cross sections for deuterium 
production via a + a — > D + • • • are unavailable. To 
make an estimate of deuteriu m production via a ct reac- 
tions we use exponential fit of lMercer et alJ 1)200111 to the 
cross sections for a + a — > d + 6 Li, the channel with the 
lowest threshold, and ignore other possible final states 
(Ad, cid 4 He). Also, the only existing cross sections for 
ct+C and a+O production of deuterium are given for the 
single energy E„ = 58MeV, or e n = 14.5 MeV/nucleon 
ijBertrand. Peelle. fc Kalbach-Clind 11974) . In this case 
we approximated the cross sections with: 

-2-KTj 



a e 




e > eth 
e < eth 



(11) 



where ap is normalized to e x perim e 
iJBertrand. Peelle. fc Kalbach-Clind Il97^) . 



erimental values 
and the 



Gamow factor rj 
barrier. 



jt r /hc(3 accounts for the Coulomb 



2.2. 



Results 

Before we present our final results let us first make an- 
alytic estimates for couple of reactions for this model. In 
order to perform analytic estimates some approximations 
will be needed. To analytically solve thick-target case we 
approximate energy losses via: 

dR „ „ , A m e (m) 
— = 0.14 — — - e 

de 2t:Z 2 z m v e 4 



(12) 
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Fig. 1. — Total yields of neutrons, deuterium and Li produced via spallation in the thick-target model, in particles produced per flare 
proton above 30 MeV. For clarity, we show only the reaction that dominates production of each element over most of the range of s. The 
points at s = 4 represent our analytic estimates for this choice of spectral index; n=circle, d=triangle, 6 Li=square. Abundances used is 
J/He = 0.1. Threshold energies and flat cross-sections were approximated by e t h ~ 25 MeV and ap a Ri 40 mb for deuterium production, 
e t h as 22 MeV and <r™ a Ri 60 for neutron production, e t ^ Ri 11 MeV and ~ 20 mb for 6 Li. 



Let us also assume that a cross-section for a particular 
reaction is flat and that we can neglect relativistic cor- 
rections which is a valid assumption for low projectile 
energies. Therefore, equation (6) now has the following 
form: 



= 5.7x10" 



(s - 3)Z 2 



(lOmb) ( 



eth 



10 MeV 



eth 



(13) 

In Figure ^ we present our analytic estimates along 
with numerical results for deuterium production by pa 
, 6 Li production by aa and neutron production by pa re- 
actions. The analytic result is calculated for the case of 
s = 4 which is the best choice given the assumptions 
that were made, while numerical results are functions of 
spectral index s. Note that our results also include reac- 
tions where target and projectile species are interchanged 
(initial-state particle species can also serve as targets, 
since we assume that flares have the same composition 
as the ambient medium) . Although we made some rough 
estimates we see that they are in good agreement with 
our numerical predictions, which gives us confidence in 
our code. 

Our numerical results for production of neutrons, deu- 
terium and 6 Li are shown in Figure Results are rep- 
resented in the form of ratios of total numbers of parti- 
cles produced in spallation as a function of the spectral 



index s. The purpose of this was to see how neutron 
production compares to production of deuterium, since 
the key idea of ML paper was that a significant amount 
of neutrons can be produced via spallation which can 
then undergo radiative capture to make a non-negligible 
amount of deuterium. We see that indeed, n/d > 1, so 
that neutron production does in fact dominate direct d 
production for all spectral indices s. More specifically, 
n/d ~ 1 around s = 4, the best-fit value for solar flares, 
and can be as high as n/p ~ 200 for an extreme value 
of s = 7. These calculations thus confirm the suggestion 
of ML that in fact n production is significant compared 
to deuteron production, and thus radiative capture of 
the neutrons offers an important channel for deuterium 
synthesis which had been neglected by ELS. 

With the neutron production in hand, we can now 
update the ELS argument to address the ML loophole. 
The key point here is that spallation production of d 
and n is also inevitably accompanied by production of 
other light elements, notably lithium. In particular, 
since 6 Li is uniquely produced by spallative pr ocesses 
(IFields fc 01iveill99a IVangioni-Flam et alJll999D . it of- 
fers the strongest constraint, as follows. Assume that all 
neutrons made in flares undergo radiative capture onto 
protons, rather than suffering decay or non-radiative cap- 
ture on 3 He. We also assume that all 6 Li in the ISM 



5 




Fig. 2. — The ratios of total numbers of secondary particles produced via spallation in the thick-target case. All ratios are shown as 
functions of the spectral index s. 



comes from flares, which maximizes the flare contribu- 
tion. Then by using our (n + d)/ 6 Li ratio for the en- 
tire range of spectral indices, 20 < (n + d) / 6 Li < 1000, 

combined with the solar 6 Li/H = 1.5 x 10~ 10 , we in- 
fer a flare-produced deuterium abundance in the range 
3 x 10~ 9 < D/H < 1.5 x 10~ 7 , much smaller than ob- 
served deuterium abundance. This limit would be further 
strengthened if we also note that a large fraction (> 50%) 
of neutrons escape the Sun, or are captured onto 3 He 
rather than on protons (Ramatv 1986). If f np < 0.5 is the 
fraction of n which do capture on protons, thus the cor- 
rect ratio for D/ 6 Li is (f np n+d)/ 6 Li = (f np n/d+l)d/ 6 Li, 
which lowers the total deuterium production by an addi- 
tional ~ 25% for n/d ~ 1 and f np = 1/2. 

We can also go the other way. Assume that just 10% 
of observed ISM deuterium abundance comes from flares. 
Then from our (n + d)/ G Li ratio for all spectral indices 
we get that 6 Li/H abundance is in the range 1.5 x 10~ 9 < 
6 Li/H < 7.5 x 10~ 8 , which is between 1 and 2.5 orders 
of magnitude larger than the observed solar abundance 
of 6 Li (where overproduction of about 300 corresponds 
to the spectral index most favored for solar flares). Thus 
we conclude that then we can rule out the ML loophole, 
if flare spallation products escape with equal probability. 
This result updates and reaffirms the same conclusion by 
ELS. 

However, many flare stars will have fewer 
metals, and thus a different composition. For 
example, disk G and K dwarfs are known 



( Ko toneva. Flvnn Chiappini. &: Matteuccil l2002t 
Wvse & Gilmore 1995) to have a mean metallicity 
that is subsolar by about a factor of 2. What effect will 
this have on the production ratios and ultimately the 
revised ELS constraints? To get a sense for this, we 
follow ELS and repeat our calculation for the extreme 
limit of a primordial composition-i.e., where the flare 
and the ambient medium contain only 4 He (with mass 
fraction 24.8%) and the balance is H. The difference 
here is that we exclude reactions involving CNO, but 
since Li can be made by aa fusion, the 6 Li constraint 
remains. In this case the (n + d)/ 6 Li ratio is in the 
range of 10 < (n + d)/ 6 Li < 30000, as shown in Figure 
thus placing the upper limit (for very low spectral 
indices) on the expected deuterium abundance due to 
flares about an order of magnitude higher. Therefore 
we see that even in primordial environment the ELS 
constraint holds, especially since thick-target production 
ratios that we are discussing change only slightly around 
most favorable spectral index of s « 4. 

Now, although our fiducial calculation is based on 
the solar abundances, the Sun itself shows a differ- 
ent and more enriched composition in flare projectiles 
and possibly targets. In solar flares, 4 He/H can be as 
much as 5 times higher than in mean solar matter- 
in which case protons and a particles can be within 
a factor of two of each other by number; abundances 
of other h eavy elements can also be similarly enhance d 
in flares ( Mandzh avidze. Ramatv. fc Kozlovskvl [l999). 
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Fig. 3. — As in Figure |21 but with projectile and target abundances as follows: (a) primordial, i.e., He/H = 0.1, no metals; and (b) 
metal-enriched, as is seen in solar flares. Here we take y a = 0.5, yo = 4.24 X 10 -3 , y N = 5.60 X 10~ 4 and y c = 1.79 X 10 -3 . 



Therefore, one might wonder how will our results change 
in that case? Figure |2h plots our results obtained as 
before, but this time for abundances that are 5 times 
greater than solar. However, even in this case our main 
conclusions stay the same: (?) the neutron channel for 
D production can dominate over spallation production 
of D, but not enough to be a significant source of non- 
primordial D, and (ii) the ELS constraint on 6 Li/D holds, 
and indeed is the strongest in this case, since overproduc- 
tion of 6 Li goes between 1 and 3 orders of magnitude. 
Therefore we see that by using a solar composition to 
describe flare-processes, our constraints were in fact gen- 
erous compared to this case. 

Of course, ML note that it is possible that D escapes 
more readily than does Li, due to its small mass. Clearly, 
the question of transport is complex, involving competing 
effects such as convection and mass loss. These difficult 
issues in magnetohydrodynamics are beyond the scope of 
this paper, but we can at least quantify the bias needed 
to avoid our updated ELS constraint. We saw that if 
interstellar D is due to flares, then 6 Li is overproduced 
by a factor that is about between 10 and 250. Thus, 6 Li 
mass loss must be suppressed relative to D by at least 
this factor in order for the ML loophole to remain open. 
While it is difficult to rule this possibility out (or in!), 
clearly this offers a strong quantitative constraint on the 
particle loss mechanism. 

Thus, we strongly constrain the ML loophole, which 
now requires a strong bias against Li escape. But if such 
a bias could be created, is the ML loophole viable? To 
further explore that scenario, we turn to the 2.22 MeV 7- 
ray line as an additional constraint on neutron captures 
in the Galaxy. 



3. GAMMA-RAY LINE CONSTRAINTS ON FLARE 
PRODUCTION OF DEUTERIUM 

In the previous section we have seen that analysis of 
particle production through spallation reactions does not 
completely answer the question of the amount of deu- 
terium produced in stellar flares. For that reason we 
will now approach this problem from a different angle. 
The reaction of radiative capture n + p — » d + 7 which 
was proposed to be potentially significant source of deu- 
terium in ML paper, has as a result, besides deuterium, a 
2.22 MeV 7-ray line. Our idea is to predict the Galactic 
7-ray intensity of that line under the assumption that 
flares produce significant amounts of deuterium. That 
way, without going into details about mechanisms that 
can transport produced deuterium into the ISM, we can 
place an upper limit on non-primordial deuterium pro- 
duction via radiative capture. 

If we denote number of 2.22 MeV 7-ray lines produced 
by a single flare with N 1 , and use N d ot to denote number 
of deuterium produced by the same flare in both radiative 
capture and spallation processes then we can write 

1 (14) 
'7 J2.2 N n 

where N n is the number of neutrons produced by the 
same star via spallation. The factor $2.2 is the efficiency 
for ncutron-to-2.22 MeV phot on conversion, averaged 
over neutron energy spectrum l|Ramatvl 11986): this in- 
cludes stopping of the 7 rays, as well as neutron escape 
and neutron "poisoning" by capture onto 3 He. Since 
deuterium is made via radiative capture as well as in 
spallation, we have Nf* = Nf M + N r d c . Let us further 
assume that all neutrons made in spallation reactions, 
N n , go into making of deuterium by radiative capture, 
that is, assume that N d c = N n . Thus the 2.22 MeV 7-ray 



N$* 
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intensity estimated with this assumption corresponds to 
the upper limit of Galactic flare production of deuterium. 
We can now rewrite equation (14) and obtain r^ ot , the 
rate of total deuterium production by a single star, as a 
function of deuterium-to-neutron spallation production 
ratio 



jntot 
1 d 



Ix. 

h.2 



l 



spall 



(15) 



What we actually want is the Galactic deuterium pro- 
duction rate by mass, M^ a \ which is a function of total 
deuterium production rate r^ ot and total number of flare 
stars in the galaxy A/"* : 



Mf = m d T*?M* 



Ll 
h.2 



1 



spall 



mdN; 



(17) 



where A/"* is the number of flare stars in the Galaxy, and 
rrid is the deuteron mass. 

We wish to find the observable Galactic 2.22 MeV j- 
ray intensity 7 7 , in terms of the 2.22 MeV emissivity q^ al , 
the total Galactic production rate of 2.22 MeV 7 rays per 
volume per steradian 



gal = r 7 n * 

where n*(r) is the number density of flare stars, 
intensity is simply a line-of-sight integral: 



qf 1 dl 



r 7 iy los 

47T 



(18) 
The 

(19) 



where A^ os = J, n(r) dl is the column density of flare 
stars along the line of sight I. 

Therefore, we can link Galactic deuterium production 
rate by mass as a function of Galactic 2.22 MeV 7-ray 
intensity by combining eas. 1171 and 1181 



or 



M 



I, 



gal 



47T 

h.2 



1 



spall 



Ni os 



1 + (d/n 



spall 



/ 2 . 2 iv los 



(20) 



(21) 



47rm ( zA/'* 

Equation l|20l) allows us to estimate the most modest 
Galactic 2.22 MeV 7-ray intensity that we expect if 
the flare stars are producing significant amounts of deu- 
terium. In order to predict the lower intensity limit, let 
us take d/n = 1 (although from Fig. 2 we see that this ra- 
tio is always smaller then 1). We also adopt the /og = 0.1 
neutron-to-7 conversion efficiency calculated bv iRamatvl 
(1986); this value is dominated by escape (followed by 
decay), and then the neutron poisoning due to the large 
n-capture cross section for 3 He. 

The Galactic distribution n(r) of flare stars (i.e., 
M dwarfs) controls both the column density and flare 
star number which appear in eqs. (|21|l and (|20|) . 
We have verified that disk M dwarfs dominate the 
calculation, while spheroid (" stellar halo") M dwarfs 
ijGould. Flvnn. k Bahcalll9"98j) do not contribute signif- 
icant numbers or column of stars . For the disk M dwarfs, 
we adopt the ijZheng et al.ll200H) distribution 

n(r,z)=n Q exp^- r J 1 " j x (22) 
[(1 - P) exp(-HMi) +/3exp(-|z|// l2 )] 



iZheng et alJ l|200l(l use Hubble Space Telescope star 
counts to determine the initial mass function in the 
M dwarf range, and find the parameters in eq. 12 I'll 
to be: local M dwarf density no — po/(m,M dwarf) = 
5.2 x 10~ 2 pc" 3 , radial scale length H — 2.75 kpc, solar 
Galactocentric distance Rq = 8 kpc, and vertical scale 
parameters hi — 156 pc, hi — 439 pc, and (3 = 0.381. 
With these, we find the number of Galactic M dwarfs to 
be A/"* = 2.4 x 10 10 , and a column density towards the 
Galactic center of Ni os = 5.1 x 10 3 pc" 2 . 

The deuterium production rate that we would consider 
to be significant is one that is comparable to the Galactic 
destruction rate of deuterium due to its cycling through 
stars ("astration"), which burn D during the pre-main 
sequence phase. The present rate at which deuterium is 
lost is 



M d , 



= -X d £ = -2.5 x lO" 5 M /yr 



(23) 



where Xd = 2.1 x 10" 5 is the present ISM mass frac- 
tion of deuterium, and £ is the current rate at which 
mass is ejected from dying stars. This is roughly 
given by £ = Rip, where R ~ 0.3 is a conserva- 
tive e stimate of the "return fraction" of mass from 
stars (Pa gellll997l e.g.), and we adopt a current Galac- 
tic star formation rate ij) ~ 4 /yr follo wing, e.g., 
iFields. Olive. Casse. k Vangioni-Flaml 1)200 lj) . 

We then can place a lower limit on the 2.22 MeV Galac- 
tic 7-ray intensity that corresponds to flare production 
of deuterium equal to the amount of deuterium that is 
destroyed per year in the Galaxy: 

J 7 = 4 x 10~ 2 cm" 2 sec" 1 sr" 1 (24) 

We can compare this predict ion with 2.22 MeV al l- 
sky map done by COMPTEL l)McConnell et all Il997|) . 
This map contains a single possible point source but no 
diffuse emission is found, down to a level estimated (M. 
McConnell, private communication) to be 



I° hs < 5 x 10" 4 cm" 2 



sec 1 sr 1 



(25) 



This limit is about 1% of what we would expect if a 
significant amount of deuterium originates in flares (eq. 
124(1 . Or to go other way around, we can use eq. I|21|l to 
translate the COMPTEL upper limit into an upper limit 
for deuterium production by flares: 

Mf 1 < 3.0 x 10" 7 M Q /yr ~ 1 x 10- 2 |M d , astrate | (26) 

Thus we see that the maximum possible D production 
rate is at most 1% of the destruction rate, so that the 
net effect is that D is indeed destroyed and thus decreases 
with time. If we were only to look at the n channel of 
deuterium production, as it was proposed by ML, we 
would get even lower production rate (by about a factor 
of two). Therefore on the basis of COMPTEL observa- 
tions and our numerical analysis of production rates we 
can conclude that stellar flares are not significant sources 
of non-primordial deuterium on the Galactic scale, and 
thus do not spoil the monotonic decline of D with time. 

Although we have estimated the diffuse 2.22 MeV ra- 
diation from flare stars, our calculation in fact is gen- 
erally applicable to any Galactic sources of deuterium 
via neutron radiative capture. The 7-ray constraints are 
thus stronger than the abundance-based constraints of 
the previous section. Furthermore, further limits (or ob- 
servations) at 2.22 MeV will constrain any source of D 



production from neutrons. We encourage INTEGRAL 
observations to be made to tighten the constraints we 
present here. 

While 7-ray observations rule out flare D production 
at a level which would overwhelm destruction, the fact 
remains that neutron production is a significant and in- 
evitable result of flare activity, as emphasized by ML. 
Thus, a 7-ray signature of this process must exist at 
some level. We can make a crude prediction of diffuse 
2.22 MeV 7-ray intensity that INTEGRAL might ob- 
serve. The rate at which 2.22 MeV 7-ray are produced 
in a single star is given by 

r 7 = /2.2-rT-kflare (27) 

where Lfiaro is the time-averaged luminosity in the flare 
state, and 

(e) = Y^Je Ni(e) de = J ek p e- s de £ Vl A t (28) 

measures the flare energy going into accelerated particles. 
Using the spectral index s = 4 favored by solar flares, 
and integrating from the lowest threshold energy among 
reactions involved in neutron production (note that be- 
cause of assumed power law spectrum of flare particles it 
is not possible to extrapolate to zero energy) we find that 
the number of neutrons produced per flare unit energy is 
Qn/( e ) = 0.015 atoms erg -1 . 

For the total number of neutrons produced per pro- 
ton (above 30 MeV), we used our thick-target s = 4 
numerical result Q n — 5.53 x 10 -3 . Then by tak- 
ing a bolometric flare lumi nosity of dM5e star to be 
L flarc = 1.9 x 10 29 erg/sec HColeman fc Worded Il976l 
and neutron-to-photon efficiency factor to be /2.2 = 0.1 
it follows from equation (19) and (20) that 

J 7 - 1.2 x 1CT 8 cm -2 s -1 sr -1 (29) 

In this cal culation we too k V* — 10 12 pc 3 and n* = 
0.86 pc~ 3 l|Weistrodll972|) . and we took the a line of 
sight length I = 20 kpc towards the Galactic center. Fi- 
nally, we note that using the same estimate of flare ac- 
tivity, (i.e., scaling to the flare star luminosity), we can 
also arrive at an estimate of the flare D production rate, 
namely M ( f are ~ 2.5 x 10~ 10 M Q /yr, which is completely 
negligible compared to the D astration rate. 

The emission in eq. (|29f) is low, and in particular is too 
dim to be observed by INTEGRAL. On the other hand, 
our estimate is crude, and while detection of flare stars 
would be a surprise, it would also provide unique new 
information about the particle content and energetics in 
stellar flares. Also, as we have noted, simply tightening 
the limits on diffuse 2.22 MeV radiation will strengthen 
the case that D production is negligible, and can go fur- 
ther to constrain (or pro be) flares as a contri butor to the 
D variation in the ISM ijHoopes et alJl2003j) . For these 
reasons, we encourage INTEGRAL observations at 2.22 
MeV, particularly towards the Galactic center where the 
emission should be the strongest. 2 

2 We note that COMPTEL IMcConnell et alJ [1991 
ISchonfelder et all 120001) found a single point source at 2.22 
MeV, with flux </> 7 = (3.3 ± 0.9) X 10" 5 cm" 2 s -1 . A flare star 
with luminosity r 7 (eq. 1271 at distance r has a time-averaged flux 



4. CONSTRAINTS ON LOCALIZED DEUTERIUM 
PRODUCTION 

It is important to note, however, that this conclusion 
applies to the Galaxy as a whole. While the global D pro- 
duction rate due to flares is small, it is a separate ques- 
tion whether D production could be sufficient to create 
D variations on smaller scales. 

To constrain local enrichment, we rely on the Galactic 
average deuterium production rate M^, which we have 
constrained in the previous subsection on the basis of 
gamma-ray observations, and estimated using flare en- 
ergetics. With the global production rate in hand, we 
can find the average D production rate for one star. The 
limit from eq. gives ra d = M d /N* ~ 10~ 17 M /yr. 
Using this, let us find the needed D production to pol- 
lute a parcel of gas at a level of 5D/D ~ 5Md/Md ~ 1, 
i.e., to product a factor of 2 variation in D/H over the 
current ISM level. To pollute a gaseous region of total 
mass M g given the entire age of the universe requires 
that A^m dwarf ~ 70M S , i.e., a gas fraction of no more 
than 1/70, that is, the M dwarfs must outweigh the am- 
bient gas by a large factor. On the other hand, a clump 
of M dwarfs in the hot and diffuse ISM could achieve 
this level. If we instead adopt the deuterium production 
rate M d ~ 10~ 10 Af©/yr, based on energetics, we find 
a needed gas fraction < 10~ 5 . This becomes harder to 
achieve, but is difficult to completely exclude. 

5. DISCUSSION AND CONCLUSIONS 

For more than a quarter century, the ELS argument 
that all deuterium is primordial has played a central role 
in cosmology. The importance of this argument demands 
that its assumptions be carefully checked. In this con- 
text, Mullan & Linsky's discovery that D production by 
neutron capture in flares can evade the ELS constraints is 
thus very important and demands further investigation. 

In this paper, we have examined and constrained the 
ML scenario in two ways. First, we have made a detailed 
calculations of spallation yields in flares. These calcula- 
tions show that flares have 1 <n/d< 200, confirming the 
ML suggestion that neutron production is at least as sig- 
nificant as direct deuteron production, and considerably 
more for some spectral indices (albeit ones atypical of so- 
lar flares). We have thus updated the ELS constraints on 
flares, now taking into account neutron production (all 
which we assume leads to deuterium formation by radia- 
tive capture). Specifically, we consider 6 Li production 
which accompanies n and d synthesis. The ELS argu- 
ment is quantitatively changed, because the 6 Li/(n + d) 
ratio is lower than the former 6 Li/c£ ratio, but the qual- 
itative conclusion remains, that the solar 6 Li abundance 
forbids a significant d component from flares. The only 
way that flare nucleosynthesis could avoid this bound is 
for d to escape the star preferentially with respect to 6 Li, 
the needed bias being a factor around 1000 in escape ef- 
ficiency. 

(4>-y) — 3 X 10 _12 (r/l pc)~ 2 cm" 2 s _1 . Thus, the point source 
is too bright to be even a nearby flare star unless COMPTEL 
happened to see a brief, intense burst, from a nearby star with 
flare luminosity much larger than the average dMe. INTEGRAL 
observations of this source can test the flare star hypothesis by 
looking for 2.22 MeV time variability in coincidence with, e.g., Ha 
and other line emission, as well as X-ray and UV signatures. 
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Second, we considered the implications of the 2.22 MeV 
7-ray line produced in the radiative capture. If D is 
currently produced in the Galaxy via this mechanism 
(either in flare stars or elsewhere) then there inevitably 
is an observable 2.22 MeV line signature. In the case of 
flare stars, the emission would be unresolved and would 
lead to a diffuse intensity tracing the Galactic plane. The 
observed limits on 2.22 MeV emission from COMPTEL 
immediately translate into a constraint on D production 
by radiative capture, and rule out this mechanism as a 
significant source of D at the global Galactic scale. 

It is more difficult to constrain flare production of D as 
the source of the possible scat ter in local D abund ances as 
reported recently by FUSE ijHoopes et alJl2003|) . How- 
ever, arguments based on energetics and on gamma-ray 
fluxes do demand that the mean D production rate per 
star is small. This means either that (I ) to produce sig- 
nificant fluctuations requires a high local concentration 
of M dwarfs relative to the gas they pollute, or (2) that 
flare D production is dominated by a very small fraction 
of M dwarfs, which dominate the global mean produc- 
tion and can thus create local anomalies. In either case, 
the nearest local sources may be detectable through their 
gamma-ray lines. 

Another constraint on the ML scenario comes from ob- 
servatio ns of deuterated molec ules towards the Galactic 
center. iLubowich et"aD pTOM detect DCN in a molec- 
ular cloud 10 pc away from the Galactic center, and on 
the basis of a molecular chemistry model estimate that 



D/H = (1.7±0.3) x 1CT 6 in this cloud. Unfortunately, the 
fractionation corrections here are large and thus the re- 
sults are somewhat mo del- dependent. With this caveat, 
the important point here is that D /H is found to be lower 
towards the Galactic center, suggesting that D indeed has 
a positive Galactocentric gradient, which argues against 
any significant stellar source of D, including flares. 

Finally, we note that while flare star radiative capture 
synthesis of D is insufficient to alter the conclusions of 
ELS, it is certain that the process does occur at some 
level-the flares exist, and must produce neutrons. Thus, 
the 2.22 MeV signature of this process must exist at some 
level. We have made a simple estimate of the surface 
brightness towards the Galactic center. We find this to 
be below the sensitivity of INTEGRAL, but given the 
crudeness of our estimate, it is worth investigation, as 
limits (or detection!) of this line has immediate impli- 
cations for stellar flares, neutron capture processes, and 
deuterium evolution in our Galaxy. 
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